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Random Matrix Theory (RMT) has elaborated successful predictions for Dirac spectra in field 
theoretical models. However, a generic assumption by RMT has been a non-vanishing chiral 
condensate E in the chiral limit. Here we consider the 2-flavour Schwinger model, where this 
assumption does not hold. We simulated this model with dynamical overlap hypercube fermions, 
and entered terra incognita by analysing this Dirac spectrum. The usual RMT prediction for the 
unfolded level spacing distribution in a unitary ensemble is confirmed to a high precision. The 
microscopic spectrum does not perform a Banks-Casher plateau. Instead the obvious expectation 

1 /3 

is a density of the lowest eigenvalue X\ which increases °c A, . That would correspond to a 
scale-invariant parameter °< AV 3 / 4 , which is, however, incompatible with our data. Instead we 
observe to high precision a scale-invariant parameter z « AV 5 / 8 . This surprising result implies a 
microscopic spectral density «Aj , which still remains to be understood in the light of RMT. 
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1. Dirac spectra and Random Matrix Theory 

Random Matrix Theory (RMT) has been applied extensively to predict the microscopic spec- 
tral densities of Dirac operators in fermionic quantum field theories. These predictions were well 
confirmed, first with staggered fermions restricted to the sector of topological charge v = [j|], and 
later with Ginsparg-Wilson fermions also for charged topological sectors |§, [|]. These applications 
referred to the case of a finite chiral condensate £ = —(yip) in the chiral limit of fermion mass 
m —>■ 0, which may occur spontaneously (as in QCD), or due to an anomaly (as in the 1 -flavour 
Schwinger model). In this case the microscopic spectrum displays a plateau near zero; its value is 
directly related to £ by the Banks-Casher formula (in finite volume the plateau is slightly shifted 
away from zero). In the £-regime RMT predicts in addition a wiggle structure superimposed on 
this plateau, which does in fact agree with numerical data. Matching them to the RMT curves 
determines £ as the only free parameter — this is a neat way to evaluate £. 

There are also models with £ = 0, but that situation is hardly explored by RMT. It occurs 
for instance in systems of fermions interacting through Yang-Mills gauge theory above the critical 
temperature of the chiral phase transition. There are numerical studies and conjectures about it 
[Q, ||], but the features of such spectra remain controversial. 

In addition there are models with £(m — > 0) = even at zero temperature. This is the case for 
the 2-flavour Schwinger model, which we are going to discuss here. We simulated this model with 
dynamical Ginsparg-Wilson fermions [Q]. Here we present our observations on its Dirac spectrum, 
which has been completely unexplored so far. 

First we review the model and the analytical predictions for £. Next we sketch our lattice 
formulation and its simulation. Then we present the unfolded level spacing distribution, which 
agrees accurately with the generic RMT prediction for the unitary ensemble. As the main subject 
of this report, we then focus on the probability density of the leading non-zero Dirac eigenvalue 
Ai, where we reveal a surprising result. The scale-invariant parameter — a rescaled eigenvalue in 
finite volume — does not agree with the obvious conjecture, which is based on the critical exponent 
8 derived in the literature. The microscopic spectrum does increase without a plateau, but its slope 
follows an unexpected power law. Finally we also consider the bulk eigenvalues. 



2. The chiral condensate in the Schwinger model 

The Schwinger model corresponds to QED in d = 2, given by the Lagrangian 

JSf(*,¥,A M )=*(*) [y^+gA^) + m]^{x)+ l -F llv {x)F llv {x) . (2.1) 

For Nf degenerated fermion flavours of mass m<g,a bosonised form of the Schwinger model 
leads to the prediction [0], 

SW« m '/ 5 , S = ^T- (2.2) 
Nf-l 

• In the quenched case (formally Nf = 0), the divergence of £(m — ► 0) agrees with simulation 
results |g|. 

• For Nf = 1 — the original version of the Schwinger model — one obtains a finite value 
£(m — > 0) = e y / (2k 3 / 2 ) due to the axial anomaly ( y is Euler's constant) 
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• The explicit prediction for Nf = 2 reads Z(m) ~ 0.38m 1//3 [10]. 

A numerical study with Domain Wall Fermions on a 16 2 x 6 lattice measured £ in the range 
m = 0.1 . . . 0.3, and a fit in this regime suggested £(m) <=c m o.388(68) [jfjlj 



3. Lattice formulation and simulation of the 2-flavour Schwinger model 

In our study, we use a lattice formulation with compact link variables U^ yX £ U(l) and the 
plaquette gauge action. For the fermions we apply the overlap hypercube Dirac operator, 



D, 



ovHF 



i-=)a&F+" 



D 



(o) 

ovHF 



1 + (D hf -1)/V( d hf-1)(£>hf-1) • (3.1) 



Duf(U) is a hypercube fermion operator [|l2|]: it is truncated perfect, and thus by construction 



approximately chiral [|13fl. In eq. ( p.l| ) it is inserted into the overlap formula [14]], which restores 



exact (lattice modified) chirality [|15p. The spectrum of D^ HF is located on a unit circle in the 
complex plane, with centre 1. Compared to the standard overlap fermion formulation — where the 
Wilson operator is inserted into the kernel — D ov hf has a better level of locality, it approximates 
rotation symmetry better, and it has an improved scaling behaviour |TJ, O, |TJ]. All these virtues 
are based on the similarity between the kernel and its overlap operator, D^j UF Dhf- 

In addition, that property also allows us to use a simplified form of the HMC force term, given 
by a low polynomial in Dhf- This reduces the computational effort for dynamical overlap fermions. 
The algorithm is kept exact by applying D ov hf to high precision in the Metropolis step at the end 
of each trajectory. 



In this way, we simulated this model at weak gauge coupling, j8 = \ t 



5, with two degen- 



erated fermion flavours of mass m = 0.01 . . . 0.24, on L x L lattices, L = 16 ... 32 [§]. Regarding 
systematic errors, the chiral extrapolation appears safe, and lattice spacing artifacts are harmless as 
well (we always deal with smooth configurations, plaquette values ~ 0.9). Finite size effects have 
to be discussed, however. To illustrate this, we show in Fig. [I] the (theoretically predicted [[!(]]) 
correlation length in the regime of the fermion masses that we simulated. The significance of finite 




Figure 1: The correlation length t, , ranging from §(m = 0.24)~ 1.69 to §(m = 0.01)= 14.03. 

size effects also implies that the distinction between the topological sectors is important (this is 
usually characteristic for the £-regime). The latter are identified by measuring the fermion index 
JT7[]. Our HMC histories contain only few topological transitions for the light masses, hence we 
performed measurements in fixed sectors. 



3 



Chiral limit without chiral condensate 



Wolfgang Bietenholz 



4. Pattern of chiral symmetry breaking 

At least in 4d Yang-Mills theory, there are only three patterns of spontaneous chiral symmetry 
breaking, depending on the fermion representation [ |T8| ] 

SU(Nf)®SU(N f )-+SU(Nf) : unitary 

SU(2N f ) -> 0(2N f ) : orthogonal 
SU(2Nf)^Sp(2N f ) : symplectic 



(complex representation) 
(real representation) 
(pseudo-real representation) . 



RMT predicts the unfolded level spacing distribution in each of these patterns [|19p. 

For a set of configurations, conf = 1 . . .N, we proceed as follows: we numerate the eigenval- 
ues in each configuration hierarchically, A? onf , i = I...L 2 . Now we put all eigenvalues (of all N 
configurations) together and order them again hierarchically, so we attach labels k = 1 . . .NL 2 . The 
normalised difference [k(Xf°f) - k(Xf° nf )]/N is the unfolded level spacing. 
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Figure 2: The cumulative unfolded level spacing density according to RMT (for different patterns of chiral 
symmetry breaking), and from our simulation data. As in QCD, we find accurate agreement with the RMT 
curve for the unitary ensemble. A tiny deviation at L = 16 disappears as the volume is enlarged. 

We consider the eigenvalues of £>^ HF , i.e. of the Dirac operator D ov hf that we used in the sim- 
ulation, after subtracting the mass. The eigenvalues with Im A, > are mapped stereographically 
onto IR + , X\ — ► — A,/2)|. Fig. || shows that the resulting level spacing distribution is in 

excellent agreement with the RMT prediction for the unitary ensemble, as it was observed before 
in QCD [Q, |2(J. We conclude that this specific RMT formula is so generic that it is not even altered 
by the absence of a chiral condensate at m = 0. 



5. The microscopic Dirac spectrum 

In infinite volume, V —> °°, the chiral condensate is given by the Dirac spectrum as 

p(A) 



dX 



X+m 



(p : eigenvalue density) . 



(5.1) 
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Along with the prediction quoted in Section 2, £ <=c m 1 / 3 , this suggests 

p(A>0)oc A 1/3 , 



(5.2) 



in contrast to the Banks-Casher plateau that one obtains in the standard setting (with E(m — ► 0) ^ 0). 
In that case, the density for the rescaled small eigenvalues XjLV is scale-invariant (at fixed mLV) 



[21]. In our case, the very general relation (A,-) °< [Vp(X > 0)] 1 implies that the parameter [22] 

£ = A,V 3/ % ? (for small A,-) (5.3) 

is expected to adapt this role, at fixed jj.r = mV 3 ^W^ — or simply at small m. is a constant of 
dimension [mass] 1 / 2 , which is (in this context) analogous to £ in the standard setting. 

Hence we probed the corresponding finite-size scaling, but it is not confirmed. Instead our 
data are in excellent agreement with a scale-invariant parameter 



Zi = XiV 5 '^W z (W z : constant of dimension [mass] 



(5.4) 



This is illustrated in Fig. |3| for our lightest fermion mass, m = 0.01, in the sectors of topological 
charge v = and |v| = 1. We also tested the behaviour if the rescaled mass is kept m const., 
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Figure 3: Cumulative densities of A,-V 3 / 4 °< (left), and of ^V 5 ^ z,- (right), for i = 1 . . .4, at mass 
m = 0.01 and topological charge v = (above) resp. | v| = 1 (below). We see that Q strongly deviates from 
scale-invariance, whereas z, obeys this property to an impressive precision. 

as an alternative to just keeping m very small. In Fig. ^ we compare (£,■) in different volumes, 
V = L 2 , L = 16 and 32, again in the sectors | v | = and 1 , for rj const. We add the corresponding 
test with (z,i) and \i z = mV 5 ^W z « const., which displays again a superior finite size scaling. 
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Figure 4: Finite size scaling for at fir °c raV 3 / 4 « const, (left) and (zi) at pt z oc mV 5 / 8 
These plots confirm again that z\ performs much better as a scale-invariant variable. 



const, (right). 



To complete this discussion, we consider even a third scenario, where the exponent of V in the 
rescaling factor is between the two cases considered so far: now the scale-invariant variable would 
be Zj = XjV 2 / 3 Wz {Wz of dimension [mass] 1 / 3 ). The behaviour of Z\ is shown in Fig. ^ (plots 
above); as in Fig. ^] we fix m = 0.01 and consider | v| = and 1. The finite size scaling quality is 
clearly better than the one of £i, but it cannot compete with z\. This third scenario belongs to a 




Figure 5: Above: Finite size scaling for Z\ oc A^ 2 / 3 , at m = 0.01 and v = (left), or |v| = 1 (right). 
Regarding scale-invariance, Z\ performs better than £i, but not as good as z\. Below: Eigenvalue histogram 
for m — 0.01, v = compared to the spectral density pAiry in eq. (5.5), which RMT predicts in this case. 
None of these three plots does convincingly support this scenario, in contrast to the compelling evidence that 
we found for that finite-size scaling of the variable z, = AjV 3 ' 5 W Z ■ 

theoretically explored universality class: it corresponds to p(A > 0) oc A 1 / 2 , which is the spectral 
density obtained by RMT in the Gaussian approximation. There is a detailed prediction for the 
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spectral density in terms of Airy functions Ai [Q], 

p A iry(Z)°<Z[Ai(-Z)] 2 +[Ai / (-Z)] 2 {~Vz/n atZ»l). (5.5) 

Fig. H (below) compares this function to the histogram that we obtained in various volumes at 
m = 0.01 and v = 0. Our data display a far more marked wiggle structure, hence this agreement is 
not convincing, and we stay with Zi = XjV 5 ^W z as the clearly preferred scale-invariant variable — 
even though no theoretical prediction for the detailed structure of p (z) has been worked out so far. 
The eiTor on the exponent 5/8 will be estimated later, see last entry in Ref. [gj. 

6. Higher eigenvalues 

At last we take a look at Aio as one of the bulk eigenvalues, and we find a optimal finite size 
scaling for X\qL , see Fig. |6] (left). The plot on the right shows that this factor works well also for 
the rescaled full cumulative density (including all eigenvalues up to the considered value). Based 
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Figure 6: For Aio, a bulk eigenvalue, the scaling factor is shifted to L 115 . The rescaled full spectral cumula- 
tive densities in different volumes (for m = 0.01, V = 0) agree well, and turn into bulk behaviour p (A) <*; A, 
(resp. Pcumuiative(A) °= A 2 ), which is expected in d = 2. 

on the fact that the spectral cutoff Amax = 2 is fixed in any volume, it is now tempting to speculate 
that the volume factor for a good finite size scaling gradually decreases from V 3 / 5 . . . V°. However, 
considering eigenvalues above the regime shown in Fig. ^, but below the cutoff regime, we could 
not find any consistent scaling factor — and indeed there is no need for it to exist. 

7. Conclusions 

We presented a pioneering numerical study of a microscopic Dirac spectrum near a chiral 
limit with L(m — > 0) = at zero temperature. In particular we analysed spectral data of the Nf = 2 
Schwinger model, obtained from simulations with dynamical chiral fermions. 

The unfolded level spacing density follows the RMT formula for the unitary ensemble. 

Refs. |0, |T§] predict £( m) « m 1 / 3 , which suggests a microscopic density p(A >0) °< A 1 / 3 , and 
the scale-invariant variable £ oc AV 3 / 4 (generally p oc A a suggests a scale-invariant oc AV 1 ^ 1+a )). 
However, this conjecture does not agree with our data. (Note that its derivation may be invalidated 



by inserting a spectral density with explicit mass-dependence, p(A,m), in eq. (5.1 ).) 

An alternative scenario with p(A >0) °< A 1 / 2 is favoured compared to the initial guess. It has 
a known theoretical background, but the data do not strongly support it either. 
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Our data strongly favour z x AV 5//8 as the scale-invariant variable, and therefore p(A >0) oc 
A 3 / 5 . This determination is much more reliable than a direct fit to the measured p(A). Interestingly, 
the first work in Ref. [ |l0| ] specifies that £ oc m 1 / 3 is expected for £ = y/lmL 3 ! 2 / '(/37T) 1 / 4 » 1, 
whereas £ -C 1 <C 2L/ \J implies £ oc m L. For m = 0.01 we are in an intermediate regime, 
£ = 0.5 ... 1 .3 (and ILj \J ttjS = 8.1... 16.2), which renders our exponent in £ oc m 3 / 5 plausible. 

For a precise theoretical test, we hope for RMT formulae to be worked out for this setting, so 
they can be confronted with our results; this is not straightforward, but it may be feasible [22]. 
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